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We theoretically consider charge transport through two quantum dots coupled in series. The corre- 
sponding full counting statistics for noninteracting electrons is investigated in the limits of sequential 
and coherent tunneling by means of a master equation approach and a density matrix formalism, re- 
spectively. We clearly demonstrate the effect of quantum coherence on the zero-frequency cumulants 
of the transport process, focusing on noise and skewness. Moreover, we establish the continuous 
transition from the coherent to the incoherent tunneling limit in all cumulants of the transport 
process and compare this with decoherence described by a dephasing voltage probe model. 
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Introduction. The analysis of current fluctuations in 
mesoscopic conductors provides detailed insight into the 
nature of charge transfer 0, . The complete informa- 
tion is available by studying the full counting statistics 
(FCS), i.e. by the knowledge of all cumulants of the dis- 
tribution of the number of transferred charges 2,3]. As a 
crucial achievement, the measurement of the third-order 
cumulant of transport through a single tunnel junction 
was recently reported p|. To what extent one can ex- 
tract informations from current fluctuations about quan- 
tum coherence and decoherence is the subject of intense 
theoretical investigations: e.g. dephasing in mesoscopic 
cavities and Aharonov-Bohm rings J5j and decoherence 
in a Mach-Zehnder interferometer [||. 

Quantum dots (QDs) constitute a representative sys- 
tem for mesoscopic conductors. Recently, the real-time 
tunneling of single electrons could be observed in QDs 
providing an important step towards an experimen- 
tal observation of the FCS. For single dots the FCS is 
known to display no effects of quantum coherence H, . 
In contrast, in serially-coupled double QDs the su- 
perposition between states from both dots causes promi- 
nent coherent effects. Noise properties have been studied 
theoretically both in the low 1 1 11 ] and finite frequency 
range f° r these structures but no FCS studies are 

available yet. Experimentally, the low- frequency noise 
has been investigated very recently in related double-well 
junctions [l4| . 

In this Letter we show that detailed information about 
quantum coherence in double QD systems can be ex- 
tracted from the zero-frequency current fluctuations. For 
this purpose we elaborate on the FCS in the limits of co- 
herent and incoherent transport through the QD system 
by means of a density matrix (DM) and master equa- 
tion (ME) description. We demonstrate a smooth transi- 
tion between these approaches by decoherence originat- 
ing from coupling the QDs to a charge detector. The 
results are compared to a scattering approach, where 
decoherence is introduced via phenomenological voltage 
probes. 

Model. The central quantity in the FCS is P(N,t ), 
the distribution function of the number TV of transferred 
charges in the time interval to. The associated cumulant 



generating function (CGF) F(x) is 

exp [-F(x)\ = J2 P ( N > *o) exp [iN X ] 



(1) 
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Here we consider the zero frequency limit, i.e. to much 
longer than the time for tunneling through the sys- 
tem. From the CGF we can obtain the cumulants 
Cfc = — (— id x ) k F(x)\ x =o which are related to e.g. the 
average current (7) = eC\/to and to the zero- frequency 
noise S = 2e 2 C2/to- The Fano factor is defined as CijC\. 
The skewness of the distribution of transferred charges is 
given by the third-order cumulant C3. 

The setup of the coupled QD system is shown as the 
inset of Fig. ^ QD1 is connected to the emitter with a 
tunneling rate T e and QD2 to the collector contact with 
rate T c . Mutually they are coupled by the tunnel matrix 
element fl. One level in each dot, at energies e\ and £2 re- 
spectively, is assumed. We consider zero temperature and 
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FIG. 1: (Color online) Current statistics for O/r = 0.5 and 
for various dephasing rates T v /T —0, 5, 20; dashed lines: 
master equation (ME) approach, solid lines: density matrix 
(DM) formalism; on-resonance Ae = 0, symmetric contact 

coupling: r = r c = r c . r = (2ra 2 )/[4^ 2 + r(r + r v )]. 

Inset: Setup of the coupled QD system with (e)mitter and 
(c)ollector contact and mutual coupling Q. 
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work in the limit of large bias applied between the col- 
lector and emitter, with the broadened energy levels well 
inside the bias window. To compare DM/ME- and scat- 
tering approaches we consider noninteracting electrons 
(spin degrees of freedom decouple, we give all results for 
a single spin direction) throughout this Letter. We note, 
however, that strong Coulomb blockade can be treated 
within the DM/ME- approaches along the same lines. 

Coherent tunneling. The FCS for coherent tunneling 
through coupled QDs can be obtained from the approach 
developed by Gurvitz and coworkers in a series of papers 
[Til Il5l| (for related work see e.g. Ref. 16]). Starting 
from the time dependent Schrodinger equation one de- 
rives a modified Liouville equation, a system of coupled 
first order differential equations for DM elements p^g(to) 
at a given number N of electrons transferred through 
the QD system at time to. Here a,j3 G {a, 6, c, d}, where 
a, 6, c and d denote the Fock-states 1 00) , 1 10) , 1 01 > , 1 1 1 ) of 
the system, i.e., no electrons, one electron in the first 
dot, one in the second dot, and one in each dot, respec- 
tively. The probability distribution is then directly given 
by P(N,t ) = p£(t ) + pUh) + p£(*o) + P&(«o). The 
FCS is formally obtained by first Fourier transforming 
the DM elements as p a p{xM) = J2n Pa/3( t o)e' lNx ■ This 
gives the Fourier transformed equation p = C c (x)Pi with 
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(2) 



and p = (paa, Pbb, Pec, Pdd,'Re[pbc},l m [Pbc}) , r = (r c + 
r c )/2, Ae = £l -e 2 . 

Note that the counting field \ enters the matrix ele- 
ments in PJl. where an electron jumps from QD2 into 
the collector contact. The CGF is then obtained as the 
eigenvalue of C c which goes to zero for x = 0, as required 
by probability conservation [see Eq. QJ] 



F C (X) = I 



2r 



pi 



16r 2 r2 2 (e l 



1/2' 



(3) 
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with p x = 2(r 2 - 40 2 + Ae 2 ) and p 2 = T 2 + 4fi 2 
for symmetric contact coupling r c = r c = T. 

Sequential tunneling. For incoherent tunneling the 
FCS can be obtained along similar lines from a ME 
for the diagonal elements of p as p = C s p, with 
P = (Paa, Pbb, Pec, Pdd)- The coefficient matrix is 



Ux) = 



/-r c r c e l * 

T e -Z Z T c e lx 

Z -(2T+Z) 

r -r c 



(4) 



with the coupling between the single-particle states given 
by Fermi's golden rule: Z = (2\fl\ 2 /F)L(Ae, 2T) with the 
normalized Lorentzian L(x,w) = [1 + (2x/w) 2 }~ 1 [T^ |. 
The CGF corresponds to the eigenvalue of the matrix 
Q which goes to zero for x = and reads 



F S {X) 

Ql/2 



to 
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(1 + iV3)?i + (1 - iV3)q 2 + 6T + 4Z 

1/3 



- ± \/ u 2 — v 3 



(5) 



with u = 8Z 3 + 9Zr 2 (l - 3e l x) and v = 4Z 2 + 3r 2 . 

Results. The probability distributions for coherent and 
incoherent tunneling obtained from the CGFs J3J| and 
respectively, in a saddle-point approximation are 
plotted in Fig. ^ for fi/T = 0.5, where the effect of co- 
herence is most pronounced. We see that the fluctua- 
tions are smaller in the coherent limit, i.e. decoherence 
generally enhances current fluctuations. In the limits of 
small inter-dot coupling Q -C T one obtains a Poissonian 
transfer of unit elementary charges and for large coupling 
f2 3> r the FCS of a single QD is recovered In these 

limits the statistics for sequential and coherent tunneling 
are indistinguishable. 

The CGF for coherent and sequential (J5J tunneling 
yield the same expression for the average current through 
the coupled QD system [13111111: 
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FIG. 2: (Color online) Average current Ci, noise C2 in units of 
^r, Fano factor C2/C1, normalized skewness C3/C1 vs. cou- 
pling £1 for various dephasing rates T^/T —0, 5, 20; Master 
equation approach (ME): dashed lines, Density matrix for- 
malism (DM): solid lines. On-resonance: Ae = 0, symmetric 
contact coupling: r — T c —T c . 



with Ti = 2ft 2 /T. The higher order cumulants Ck with 
k >2 deviate for intermediate ft reflecting their sensitiv- 
ity to quantum coherence in the transport process. For 
r = r e = r c and Ae = we have the Fano factors 
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r 4 - 2T 2 n 2 + 8ft 4 



2e(I) 

for the coherent case and 



(r 2 + Awy 



(7) 



2e(7) 



r 4 + 2r 2 ft 2 + 8ft 4 
(r 2 + 4ft 2 ) 2 



(8) 



for the sequential, incoherent case. Clearly, coherence 
suppresses the noise [12J, [l3|. The noise and the Fano 
factors are shown in Fig. [5] (results for T v =0). The noise 
for coherent tunneling shows a local minimum at 2ft = 
r. At this coupling the normalized skewness has a local 
maximum as it can be seen in Fig.[21and a close inspection 
reveals a FCS identical to a Poissonian transfer of quarter 
elementary charges: F(x) — to^(e lx ^ 4 — 1). 

Decoherence - charge detector. In order to connect 
the limits of coherent and incoherent charge transport 
through the QD system we consider the exponential 
damping of the off-diagonal elements in the modified Li- 
ouville equation with rate T v : i.e. in the last two rows 
of the coefficient matrix T is replaced by r + Y v . 
This apparent phenomenological treatment of decoher- 
ence can be substantiated, e.g., by the introduction of a 
quantum point contact close to one of the QDs: whenever 
an electron enters the QD the transmission through the 
quantum point contact changes. This charge detection 
leads to the exponential dampin g of the off-diagonals, as 
microscopically derived in Ref. |2jj. Due to the finite 
coupling ft, it also leads to an exponential relaxation of 
the diagonal density matrix elements. Its effect on the 
FCS is presented in Fig. ^ and its effect on the current 
and noise in Fig. [21 For comparison with the sequential 
tunneling cumulants the broadening of the resonance due 
to the coupling to the quantum point contact has to be 
considered and therefore the replacement r — > V + T v in 
Z of the coefficient matrix (0} is carried out. Then, the 
currents C\ in both treatments agree for any T v (Fig. |3J) • 
The higher-order cumulants merge for T v 3> ft as shown 
for the noise C2, the Fano factor C2/C1 and for the nor- 
malized skewness C3/C1 in Fig. [21 

Decoherence - Voltage probe model. The coherent 
FCS in Eq. (j3Jl can also be obtained from the scat- 
tering formula of Levitov and coworkers |3J, F(\) = 




FIG. 3: Fano factor vs. inter-QD coupling £1 for various 
dephasing rates T v . a) elastic voltage probe, b) inelastic 
voltage probe in scattering formalism (solid curves). Dashed 
curves: master equation (ME) Fano factor for T v /T = 20. 
On-resonance: Ae = 0, symmetric coupling: r = T c = T c . 



(t /h) J deln[l + T(s)(e ix - 1)], where T(e) is the trans- 
mission probability through the QD system (see e.g. 
[Tl|). This makes it interesting to compare dephasing 
within the DM approach with dephasing in a scattering 
formalism. This is done by introducing phenomenological 
voltage probes Q coupled with strength T v = T^i = T v2 
to the QDs (see inset of Fig. |3Ji). The probes absorb 
and subsequently re-emit electrons, thereby randomizing 
their phases. Here we focus on the current and the noise, 
higher cumulants can be investigated with a modified ver- 
sion of the stochastic path- integral technique in Ref. [2l| , 
but this is beyond the scope of the present Letter. 

The scattering matrix s for the four terminal QD-probe 
system is given by 



s = 1 - iW 1 GW, G=[e-H + iWW 1 } 



Ti-l 



(9) 



H 



£1 ft 
ft £2 




The average current in lead a — e,c,ipi,ip2 is given by 



(10) 
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with Ajg 7 (e) = 8 a f}8 ai — s s af3 {e)s ai {e) and the distribution 
function f a (s) of terminal a. The zero-frequency noise 
between terminal a and (3 reads 0, |2^| 

S ^ = ^T,j deA« s (e)A Sy (e)f 7 (e)[l-Me)] (11) 

7(5 

We first consider an elastic, purely dephasing voltage 
probe [23, where the average current as well as the 
low-frequency current fluctuations into the probe is zero 
at each energy. The conservation of average current 
gives the average distribution functions f v \/ V 2- From 
the conservation of the current fluctuations one obtains 
the fluctuating part of the distribution functions Sf^i/^ 
in terms of the bare current fluctuations [lj. The to- 
tal noise is then obtained as a weighted sum of the bare 
current correlations in Eq. |H) . It is found that both 
current and noise qualitatively reproduce the DM result. 
The Fano factor is plotted in Fig. however, there is 
a quantitative difference. Since in the DM approach, the 
electrons in the dots can exchange energy with electrons 
at the quantum point contacts, the dephasing is inelastic 
and a quantitative agreement with an elastic scattering 
dephasing approach is not to be expected. 

To account for inelastic dephasing we next consider in- 
elastic voltage probes which conserve only total, energy- 
integrated current and fluctuations. Trying to mimic the 



effect of the point contacts in the DM approach, we as- 
sume the distribution functions in the probes to be con- 
stant, independent of energy in the entire bias window. 
The average current and noise are then obtained along 
the same lines as for the purely dephasing probe. We 
find that the average current coincides with the DM re- 
sult, the noise, however, again differs quantitatively but 
not qualitatively. The Fano factor is plotted in Fig. |3t>- 
We thus conclude that in double QD systems, dephasing 
in a scattering and a DM approach yield qualitatively 
similar but in general quantitatively different results. 

Conclusions. Within density matrix and master equa- 
tion approaches, we have examined the FCS for coherent 
and sequential charge transport through coupled QDs. 
While the average currents in the two cases coincide, all 
higher cumulants differ, clearly demonstrating the sen- 
sitivity of the charge transport to quantum coherence 
which generally suppresses the fluctuations. Coupling 
the QDs to a charge detector introduces decoherence, 
which results in a continuous transition from coherent 
to sequential tunneling. A scattering approach, where 
decoherence is introduced via phenomenological voltage 
probes, gives qualitatively similar results. 
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